Welcome to HySoP’s webpage!¶
HySoP (Hybrid Simulation with Particles) is a library mainly dedicated to high performance numerical simulation of fluid flows problems based on semi-lagrangian particle methods, for hybrid and distributed architectures providing multiple compute devices including CPUs and GPUs.
Source code is available at: https://gricad-gitlab.univ-grenoble-alpes.fr/particle_methods/hysop
Gallery¶
Passive scalar advection
Wake of a sphere
Sediment flows
Dissolution of carbonate rock
Playlist of several simulations
Programming elements¶
The high level functionalities and the user interface are mainly written in Python using the object oriented programming model.
This choice was made because of the large software integration benefits it can provide. Moreover, the object oriented programming model offers a flexible framework to implement scientific libraries when compared to the imperative programming model. It is also a good choice for the users as the Python language is easy to use for beginners while experienced programmers can pick it up very quickly. The numerical solvers are mostly implemented using compiled languages such as Fortran, C/C++, or OpenCL for obvious performance reasons. It is also possible to implement numerical algorithms using directly Python, which is an interpreted language, hence slower for critical code paths under heavy arithmetic or memory load. The Python language support is however the key for rapid development cycles of experimental features. It also allows to easily implement routines that compute simulation statistics during runtime, relieving most of the user post-processing efforts and enabling live simulation monitoring.
The compiled Fortran and C++ compute backends allow us to integrate a variety of external dependencies by connecting them to the main HySoP python module with interface wrappers such as F2PY or SWIG. Note that many scientific libraries already provide Python interfaces so that they can be directly used in python without needing the user to implement his own wrapper. In addition to the compiled languages, the library offers the possibility to compile generated code just-in-time during execution. This is the case for OpenCL, the language used to drive OpenCL-compatible accelerators, like GPUs, but also to translate python methods to fast machine code by using the Numba just-in-time compiler.
References¶
HySoP has been used to produce numerical results in the following references:
Jean-Baptiste Keck, Georges-Henri Cottet, Eckart Meiburg, Iraj Mortazavi, and Christophe Picard. Double-diffusive sedimentation at high Schmidt numbers: Semi-Lagrangian simulations. Physical Review Fluids, 6(2):L022301, February 2021. URL: https://hal.archives-ouvertes.fr/hal-02982080, doi:10.1103/PhysRevFluids.6.L022301.
Chloé Mimeau, Simon Marié, and Iraj Mortazavi. A comparison of semi-Lagrangian Vortex method and Lattice Boltzmann method for incompressible flows. Computers and Fluids, 224:104946, April 2021. URL: https://hal.science/hal-03214005, doi:10.1016/j.compfluid.2021.104946.
Chloé Mimeau and Iraj Mortazavi. A Review of Vortex Methods and Their Applications: From Creation to Recent Advances. Fluids, 6(2):68, 2021. URL: https://hal-cnam.archives-ouvertes.fr/hal-03677567, doi:10.3390/fluids6020068.
Jean-Matthieu Etancelin, Peter Moonen, and Philippe Poncet. Improvement of remeshed Lagrangian methods for the simulation of dissolution processes at pore-scale. Advances in Water Resources, 146:103780, October 2020. URL: https://hal.archives-ouvertes.fr/hal-02965451, doi:10.1016/j.advwatres.2020.103780.
Chloé Mimeau, Iraj Mortazavi, and Georges-Henri Cottet. Passive control of the flow around a hemisphere using porous media. European Journal of Mechanics - B/Fluids, 65:213–226, September 2017. URL: https://hal.archives-ouvertes.fr/hal-01483400, doi:10.1016/j.euromechflu.2017.03.002.
Georges-Henri Cottet, Jean-Matthieu Etancelin, Franck Pérignon, and Christophe Picard. High order Semi-Lagrangian particle methods for transport equations: numerical analysis and implementation issues. ESAIM: Mathematical Modelling and Numerical Analysis, 48(4):1029–1060, July 2014. URL: https://hal.archives-ouvertes.fr/hal-00991150, doi:10.1051/m2an/2014009.
Jean-Baptiste Keck. Numerical modelling and High Performance Computing for sediment flows. Theses, Université Grenoble Alpes, December 2019. URL: https://tel.archives-ouvertes.fr/tel-02433509.
Chloé Mimeau, Georges-Henri Cottet, and Iraj Mortazavi. Direct numerical simulations of three-dimensional flows past obstacles with a vortex penalization method. Computers and Fluids, 136:331–347, September 2016. URL: https://hal.science/hal-01855265, doi:10.1016/j.compfluid.2016.06.020.
Chloé Mimeau. Conception et mise en oeuvre de méthodes vortex hybrides-frontières immergées pour des milieux solides-fluides-poreux. Application au contrôle passif d'écoulements. Theses, Université Grenoble Alpes, July 2015. URL: https://tel.archives-ouvertes.fr/tel-01178939.
Jean-Matthieu Etancelin. Couplage de modèles, algorithmes multi-échelles et calcul hybride. Theses, Université de Grenoble, December 2014. URL: https://tel.archives-ouvertes.fr/tel-01094645.
J.-B. Lagaert, G. Balarac, and G.-H. Cottet. Hybrid spectral-particle method for the turbulent transport of a passive scalar. Journal of Computational Physics, 260:127–142, 2014. URL: https://www.sciencedirect.com/science/article/pii/S0021999113008334, doi:https://doi.org/10.1016/j.jcp.2013.12.026.
Sarah Perez. Particle and Data-Driven Approaches for Reactive Micrometric Processes : Application to CO2 Mineral Storage with Uncertainty Quantification. PhD thesis, Université de Pau et des Pays de l'Adour, 2023. URL: https://theses.hal.science/tel-04415185 (visited on 2024-05-02).
Jean-Matthieu Etancelin, Peter Moonen, and Philippe Poncet. Improvement of remeshed Lagrangian methods for the simulation of dissolution processes at pore-scale. Advances in Water Resources, 146:103780, 2020. URL: https://hal.science/hal-02965451 (visited on 2021-03-22), doi:10.1016/j.advwatres.2020.103780.
Sarah Perez, Jean-Matthieu Etancelin, and Philippe Poncet. A semi-Lagrangian method for the direct numerical simulation of crystallization and precipitation at the pore scale. 2024. URL: https://hal.science/hal-04692528 (visited on 2024-09-19), doi:10.48550/arXiv.2409.05449.
